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ABSTRACT 

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary 
disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces 
the fragmentation of surrounding planetesimals, which generaUy inhibits further growth of the core by 
removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of 
the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations 
are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the 
fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, 
kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. 
Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's 
formation was completed when the local disk mass has already decayed to a value comparable to or 
less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or 
shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of 
envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the 
core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the 
core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based 
on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter 
can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward 
migration and larger core mass of Saturn than that of Jupiter. 

Subject headings: planets and satellites: formation — solar system: formation — planets and satellites: 
individual (Jupiter, Saturn) 



1. INTRODUCTION 

Since Jupiter resides at 5.2 AU, it is likely that Jupiter 
did not underg o significant type II migration (e.g., 
Ilda fc LinI [2001 . The mass of the disk surrounding 
Jupiter decayed to a level comparable to or less than 
Jupiter mass before Jupiter completed its formation 
and opened up a gap. Since the disk could not push 
Jupiter, its type II migration is negligible. The forma- 
tion of Saturn becomes problematic along this line. With 
an assumption that cores grow through the coUisional 
accretion of surrounding planetesimals, the accretion 
timescale for Saturn's core may be 5-10 times longer than 
that for Jupiter's core. Therefore, since Jupiter forms on 
a timescale comparable to a disk depletion timescale, a 
few Myrs, Saturn's core formed on a timescale of 10- 
30 Myrs. However, when the core formed, the disk 
mass should have been so severely depleted that the core 
cannot accrete disk gas as massive as the present enve- 
lope mass. To reconcile the inconsistency, Saturn's core 
should have formed on a timescale shorter than several 
Myrs after Jupiter's formation. Note that a scenario by 
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iWalsh et al.l (|2011[ ) wher e Jupiter first mi grates inward, 
and then outward (e.g., iMorbidelli &: Cri da 2007), also 
requires rapid formation of Saturn. 

A single planetary embryo is formed from collisional 
coagulation of planetesimals in an annulus of the disk 
along the embryo's orbit and further grows through col- 
lisions with surrounding remnant planetesimals. An 
embryo reaches the critical core mass, ~ 10 Mm, to 
start gas accretion for Saturn formation ([Mizun o 198(J 
IBodenheimer fc PollackllTOSl llkoma et alj 120001 ). How- 
ever, since massive embryos enlarge the random mo- 
tion of planetesimals, collisions between planetesimals 
are destructive. Embryos also grow by the accretion of 
fragments that result from such collisions. The eccen- 
tricities and inclinations of fragments are well damped 
by gas drag, and fragment accretion accelerates embryo 
growth. On the other hand, since such small fragments 
quickly drift inwards, the final mass of an embryo is much 
smaller than the critical core mass to start gas accretio n 
(|Kobavashi et al.l[2010l[20Tlt lOrmel fc Kobavashill2012D . 
However, Jupiter may have opened up a gap in a disk 
to truncate gas infiow to Jupiter at the present Jupiter's 
mass, which forms a radial pressure maximum near the 
edge of the gap. Radial drift o f fragments is stalle d at 
around the pressure maximum (jAdachi et al.l 119761 ). In 
addition, type I mig ration of a core is also stalled there 
(jTanaka et all [2002t [Masset et al.l [20061 ). Since a core 
grows via fragment accretion without the loss of frag- 
ments and the core itself, rapid formation of Saturn's 
core may be realized under these conditions. 

In order to investigate the rapid formation of Saturn 
core near the edge of the gap, we perform simulations. 
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which can derive accurate solution s in both hmits domi- 
nated by collisional fr agmentation pCobavashi fc Tana^ 
or coagulation (jKobavashi et al .11201 Of ). In Section 
O we investigate the radial drift of bodies near the pres- 
sure maximum. We model a disk for a simulation and 
briefly explain the method of the simulation in Section 
[21 The core growth derived by the simulations is shown 
in Section H] We discuss our findings in Section [5] and 
present our conclusion in Section [6l 

2. RADIAL DRIFT NEAR PRESSURE MAXIMA 

The velocity distinction between the gas rotational ve- 
locity Wgas and the Keplarian velocity f k is given by 
vk — Wgas = rjVK with (jAdachi et al.lll976[ ) 
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where pg^s and Pgas are, respectively, the gas density and 
pressure of disk midplane, JIk is the Keplarian angular 
velocity, and a is the distance from the sun. Equation ([T]) 
shows = at a pressure maximum. 

Gap opening by Jupiter changes the radial pro- 
file of g as surfac e de nsity, which determines 77. 
ITani gawa fc Ikoma l|2007| ) presented an analytical for- 
mula for the gas surface density Sg around the gap a^ 
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where aj is the orbital radius of Jupiter, the viscosity 
coefficient v is given by ac^fl^^ with the sound velocity 
c and a parameter a, Sg 00 is the gas surface density far 
away from Jupiter, and Mj and Mq are, respectively, the 
masses of Jupiter and the sun. We assume Sg,oo cx a~P 
and c oc a"*. Equations ([T]) and ([2]) then give ry near the 
pressure maximum formed by Jupiter as 



3Pa 



(a - ajY 



(4) 



Note that Equation ([2]) gives the density profile deter- 
mined by viscous torques. The thermal instability of 
the density profile (Rayleigh criterion) r egulates the den- 
sity profile in the vicinity of Jupiter (llda fc LinI 120041 : 
iTanigawa fc Ikomal [20071: iCrida et al.l [20061). However, 
for the purpose of finding the location of the pressure 
maximum, we can use Equation ^ as long as c/Qk ^5 
O.laj and a < 0.02. 

Bodies lose their angular momenta through gas drag 
and then drift inward, resulting in their radial transport. 
The gas drag force acting on a body with radius r and 
mass m for the relative velocity u between the body and 
gas is expressed by CD7rr^pgasW^/2 with the dimension- 
less gas drag coefficient Cd, which is determined taking 
into account Epstein and Stokes regimes as well as con- 
stant Cd for high Reynolds number. The stopping time 



® ICrida et al.l II2006I ) derived a similar formula for the density 
profile of gap. The locations of the pressure maximum (8.9 AU for 
a = 6.9x10"^; 7.2 AU for a = 10'^, where we adopt c/n^a » 0.05 
at a = aj ) arc almos t the same as the values derived from the 
formula by rTanigawa fc Ikomal II2007I '). 



normalised by is given by 

Tstop = a^l^r/u, 

where 

T = 
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We do not consider very small fragments of fstop ^ 1 
because radial drift of these fragments is negligible due 
to strong coupling with gas. The relative velocity u then 
depends on eccentricity and inclination of a body as well 
as on rj. Eccentricities and inclinations of bodies have the 
Rayleigh- typ e distributions with dispersions e* and i* , 
respectively (jida fc Makinol[l992[ ). The averaged value 
of u is estimated as {e* +i* +ri)aflK- The drift velocity Vd 
is analytically given in the limit of fstop ^ 1 or fstop 3> 1 
fAdachi et al. 1976). Averaging the solutions for bodies 
with eccentricities and inclinations in Rayleigh-type dis- 
tribution and connecting both these limits, we then have 
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Equation ([7|) covers both limits in fstop ^ 1 and fstop S> 
1. Although the drift velocity becomes very low at the 
pressure maximum because of 77 = 0, bodies with high e* 
or i* can drift due to the cubic terms of e* and i*. 

A massive embryo produced in an annulus of the disk 
stirs surrounding smaller bodies. The stirring and damp- 
ing due to gas drag determine e* and i* of bodies, de- 
pending on their radii. Since e* and i* increase with 
embryo growth, bodies drift inward from the pressure 
maximum. However, since rj is negative inside the pres- 
sure maximum, the inward drift due to e* and i* is can- 
celed out at a certain distance by the outward one caused 
by negative 77. For a = 6.9 x 10~^, p — 3/2, and q — 1/4, 
the radial distribution of Sg and rj is shown in Fig. [T][3; 
the pressure maximum is then located at 9AU. Fig. [1] 
also shows the location of Wd = as a function of the ra- 
dius of bodies. The location of Wd = for small bodies is 
close to 9 AU, because e* and i* are strongly damped by 
gas drag. Although the location of Wd = moves inward 
for a massive embryo because of its strong stirring, the 
feeding zone of the embryo, of whic h the width i s given by 
10 times Hill radii of the embryo ([Kokubo fc Ida 1998), 
is so wide that the embryo can accrete kilometer-sized 
or smaller bodies. Although larger bodies drift out of 
the feeding zone, the drift timescale is much longer than 
the timescale of collisional fragmentation among them or 
of the accretion onto the embryo. Note that even if such 
large bodies drift out of the feeding zone, small bodies re- 
sulting from collisional fragmentation of the large bodies 
drift outward and halt inside the feeding zone. Therefore, 
the assumption of Wd ~ is valid to investigate embryo 
grow at the pressure maximum. 

^ Ilnaba et al.l 11200 II ) modified terms related to rj in the analytical 
formula for fstop S> 1. We apply the modified for mula for terms 
relate d to jj but the cubic terms of e* and i* from lXdachi et al.l 

* We discuss the applied a value below. 
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Fig. 1. — The radial distribution of the surface density Eg (top) 
and ri given by Equation ^ (middle) and the location of = vs. 
the radius of bodies (bottom). In order to determine the location, 
eccentricities of the bodies are given by the gas drag damping and 
the stirring of an embryo with 0.1 (dotted lines), 1 (solid line), 10 
(dashed lines) Earth masses. The vertical lines indicate the inner 
edge of the feeding zone of the embryo residing at 9 AU. 



The location of the pressure maximum depends on the 
turbulent strength, a value. Magneto-rotational instabil- 
ity may cause turbulence in the disk; its activity depends 
on an ionization degree and the strength of magnetic 
field. Because of the low ionization degree, a magneti- 
cally decou pled midplane dead zone emerges at < 20 AU 
(|Sano et al . 20000- In the dead zone, the a value is 
expected to be as low as 10~^ -10~'^ if the net vertical 
magnetic field is we ak enough iQk uzumi fc Hirosd[20Tll: 
iGressel et alj 120121 ): the pressure maximum is then lo- 
cated at around 9 AU. For larger a f» 10"'^, the pressure 
maximum moves to about 7AU. However, the behavior 
of bodies is similar; bodies smaller than 1-10 km halt 
inside the feeding zone of an embryo formed at the pres- 
sure maximum and thereby these bodies contribute to 
the embryo growth. Therefore, the location of pressure 
maximum (or choosing a value) insignificantly affects the 
Saturn core formation. We thus fix the pressure maxi- 
mum at 9 AU in the following. 

3. MODEL 

We introduce a disk model for the initial surface mass 



density of solids Eg.o and gas Sg.o such that 
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^ Although a reduction of the vertical optical depth of the disk 
due t o dust depletion could make the dead zone small l lSano et al.l 
I200CI1 . in our model very small bodies supplied by coUisional frag- 
mentation, whose total mass is negligible as shown in Section |4l 
can produce a sufficient optical depth. 
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where Emmsn,s = 1.1 gem ^ and SMMSN,g = 89 gem ^ 
are, respectively, the solid and gas surface densities at 
9 AU in t he mi nimum-mass solar nebula (MMSN) model 
(jHayashi 1981*). We vary scaling factors Xg and Xg to find 
the conditions for Saturn core formation. Solid surface 
density evolves in the simulation but gas density is set 
to be constant with time. We consider the disks with 
Xg ~ 1 which are small enough for Jupiter to reside at 
the current position without significant type II migration 
and large enoug h to form cores for Jupiter and Saturn 
(|Ida fc Linll200a) . 

We consider a disk from 8.6 AU to 35 AU. The disk 
is divided into 14 annuli whose widths are given by 0.1 
times the distance from the annulus center to the sun. 
The innermost annulus centered at 9 AU is assumed to 
be the pressure maximum zone. As shown in Section 
[21 kilometer-sized or smaller bodies stay in the feeding 
zone of an embryo residing at the pressure maximum. 
In the MMSN disk, drift time required for larger bod- 
ies to go out of the feeding zone of an embryo with 
mass < IOM0 is longer than the formation timescale 
of Saturn's core, 1-2 x 10^ years, if the bodies are larger 
than about 30 km. For bodies with intermediate sizes, 
the drift time i s longer than their colhsion al fragmenta- 
tion timescale (|Kobavashi fc Tanakal I20T0I) . Therefore, 
we ignore radial drift of bodies in the innermost an- 
nulus = 0). The values of Sg and 77 are taken 
from those at the radial center of each annulus, given 
by Equation (P and 77 = 5.4 x 10-3(a/9AU)i/2 except 
for the innermost annulus that has 77 = 0. Each an- 
nulus has a lot of mass bins to calculate the evolution 
of the mass distribution of bodies via colhsional coag- 
ulation and fragmentation among bodies. Colhsional 
outcomes are assumed to be scaled by the dimension- 
less impact energy (j) = mim^v'^ /2{mi + m2fQn ac- 
cordin g to the simple model that 
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(|201C0 and iKobavashi et al.l (|2010( ) constructed consis- 
tently with laboratory experiments and colhsional sim- 
ulations, where 771 1 and m2 are the collider masses, v 
is the colhsional velocity between the colliders and 
is the specific energy needed to eject half of the total 
mass of the colliders. The colhsional cross sections of 
bodies depend on their random velocities determined by 
e* and i* . Since the stirring of large bodies increases 
the velocities, the mass distribution of bodies affects the 
velocity evolution. Therefore, we calculate the mass and 
velocity evolution s imultaneously using the collision rates 
(|Inaba et al.ll200ll ) and the vel ocity evolution rates due 
to mutual interac tion of bodies (lOhtsuki et al.l 2002) and 
due to gas drag (Adachi et alJll976D. With" a use of an 
upwind scheme, we calculate radial transport of bodies 
due to radial drift given by Equation ^ except for the 
innermost annulus. 

In addition, we include the enhancement of collisions 
with embryos due to their tenuous atmospheres and due 
to gas drag for small bodies with fstop ~ 1, where 
"^stop = 1 corresponds to several tens centimeter in ra- 
dius at 9A U. We apply the atmo spheric enhancement 
modeled by llnaba fc Ikomal (|2003l ) with the grain de- 
pletion factor / = 1 X 10~^ and the collision rate be- 
tween embryos and small bodies with fstop ~ 1, given by 



4 



Kobayashi, Ormel, Ida 



lOrmel fc Klah^ ()2010[ ). For the Ormel & Klahr model, 
we adopt the headwind velocity of gas to an embryo as 
{rj + e-E + iE)vK, where ce and ie are, respectively, ec- 
centricity and inclination of the embryo. 

Since small bodies are strongly coupled with gas flow, 
the model of Ormel & Klahr gives a very long timescale 
for collisions between very small bodies and embryos. In 
addition, because the collisional cascade caused by em- 
bryo growth stalls at a body size of abo ut 1-10 m, the to- 
tal m ass of smaller bodies is negligible (jKobavashi et al.l 
I20T0I ). Therefore, we neglect bodies with radius smaller 
than rmin, where rj^in — 0.01cm in our simulation. In- 
deed, embryo growth is independent of rmin if r^hi ^ 
1 cm. 



4. SATURN CORE FORMATION 

We perform simulations for core formation and growth 
starting from a monodisperse mass population of plan- 
etesimals of mass mo and radius vq with e* — 2i* = 
(2mo/3MQ)i/3. In Fig. [2 the evolution of the size 
distribution due to collisions is shown for ro = 1 km 



(mo = 4.2 X 10^^ g) in the MMSN model (xg 
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Collisional coagulation of planetesimals produces a small 
number of large bodies, planetary embryos. As embryos 
grow via the accretion of a swarm of planetesimals, mas- 
sive embryos stir planetesimals and collisional fragmenta- 
tion among planetesimals turns planetesimals into small 
fragments. Since gas drag suppresses the random mo- 
tions of very small bodies, fragmentation by collisions 
among such small bodies no longer occurs. The grinding 
down of bodies through collisions halts at ^ 10 m, where 
bodies accumulate (see Fig. [2]) . The mass distribution is 
given by three discrete components that are the largest 
bodies (planetary embryos), initial sized or slightly larger 
bodies (planetesimals), and small bodies with 1-100 m. 

The gas drag coefficient Cd is constant for kilometer- 
sized or larger bodies but it increases with decreasing 
size in the Stokes regime, given by Cd — 5.5dg/ur 
(jAdachi et al.lll976| l. Here, Ig = /g,o/Pgas is the mean 
free path of gas molecular with Ig^ = 1.7 x lO^^gcm^^. 
The fragments accumulating at 1-100 m are controlled 
by the Stokes drag. The equilibrium eccentricity e*^ in 
the S tokes regime is estimated as (jKobavashi et al.|[2010l 

MM 



(10) 



where Hm = (M/SMqY^^ is the reduced Hill radius of 
an embryo with mass M, (Pvs) = 73 is the dimensionless 
stirring rate for e* ^ and 6 = 10 is the separation of 
embryos divided by their mutual Hill radius. Since r oc w 
in the Stokes regime, e*^ is independent of u; the eccen- 
tricity is the same as that in the case rj ^ 0. When 
a body of mass m collides with a similar-sized body, 
the impact energy is given by m(e*qWk)^/4. If the en- 
ergy is much smaller than 2mQ^, the collision no longer 
produces significant fragments. The collisional cascade 
stalls around (ejJqUk)^ = ChQr, with a constant Cl ^ 1 
and small bodies accumulate around the radius satisfy- 
ing the condition, which is obtained from the condition 
and Equation pUj) as 
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Fig. 2. — Evolution of surface number and mass densities of bod- 
ies larger than radius r at 9 AU, starting from a swarm of planetes- 
imals with initial radius ro = 1 km in the MMSN model. Crosses 
indicate the initial condition and lines are size distributions at dif- 
ferent times. 
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where T is the disk midplane temperature and the 
value of Qn for a bod y of radius r = 1 m is applied 
()Benz fc Asphaud 119990 . Indeed, rf is consistent with 
the result of simulations for Ch ^ I (see Fig. [2]). In- 
terestingly, Tf is independent of 77 and Pgas because 
gas drag is determined in the Stokes regime. In ad- 
dition, Cf* (ClQd)1/Vwk « 3 X 10-3(a/9AU)i/2 is 
smaller than Hm = 1O~^(M/M0)^/'^, resulting in a very 
large accretion rate of the fragments onto an embryo 
(jlda fc Nakazawalll989l ). 

Fig. 13] shows the evolution of embryo massed at 
9AU for the MMSN disk with ro = 1km and 10 km. 
Embryos initially grow exponentially with time and the 
growth time scale is proportional to ro /Sg (jOrmel et al.l 
I2010airbl ). which is called runaway growth. Massive em- 
bryos formed by the runaway growth stir remnant plan- 
etesimals and then slower, oligarchic growth starts for 
embryos larger than ~ 10"'^ M0. Such massive embryos 
then activate collisional fragmentation of planetesimals. 
The very low e* and i* of the small resulting fragments 
damped by strong gas drag tend to accelerate the embryo 
growth. On the other hand, the gas drag also leads to 
the radial drift of t he fragments, which generall y limits 
the embryo growth (jKobavashi et al.ll201Ci |2011[ ). How- 

^" In our simulation, we define "runaway bodies" which cannot 
col lide with each other due to their large orbital separations (see 
IKo bavashi ct al. 2010). Here "embryo masses" represent the av- 
eraged mass of runaway bodies, which roughly corresponds to the 
mass of the largest body. 
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ever, near the pressure maximum due to gap opening 
by Jupiter, embryos effectively accrete such fragments 
without the loss of the fragments. In Fig. [3l the embryo 
growth starts being accelerated at ~ COlAf^ (at about 
3 X 10^ years for tq = 1 km and at about 3 x 10^ years 
for To = 10 km), because fragment production becomes 
efficient at such an embryo mass. At > Mq, the accel- 
eration slows down due to depletion of the total amount 
of surrounding bodies by the embryo accretion. Nev- 
ertheless, the embryo keeps growing beyond the "isola- 
tion mass" , which is several Earth masses for the MMSN 
model, because fragments produced in outer annuli drift 
into the region at ^ 9 AU. Eventually, an embryo reaches 
10 Earth masses in 2 x 10^ years for tq = 1 km, while the 
embryo mass is about an Earth mass in 5 x 10® years for 
To = 10 km. 
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Fig. 3. — Embryo growth at 9 AU for ro = 1 km and 10 km in the 
MMSN model (xg = Xs = 1). For reference, dotted Unes show the 
results in the case of an unperturbed disk (i.e. without Jupiter's 
gap). 

In the oligarchic growth of embryos, embryos have 
orbital separations of 10 Hill radii of themselves 
(jKokubo fc Idalll998D . Since 10 Hill radii of an embryo 
with IOM0 at 9 AU, about 2 AU, is larger than the pres- 
sure maximum region (see Fig.[T]), only an embryo larger 
than IOM0 can be formed in the pressure maximum. In 
case that the pressure maximum does not emerge (dot- 
ted lines in Fig. ^ , embryo growth is not so rapid and 
stalls at about O.IMq, because surrounding planetesi- 
mals are lost due to combination of planetesimal frag- 
mentatio n with the radial drift of r esulting fragments 
(see also lKobavashi et al.l[2010l [20TTh . Since there is no 
pressure maximum in the outer disk, embryos beyond 
10 AU cannot become larger than a body of about Mars 
mass. Meanwhile the fragmentation of planetesimals and 
the ensuing radial drift of resultant fragments remove the 
solid mass reservoir in the outer disk, which facilitates 
the formation of a large enough embryo in the pressure- 
maximum region. 

Although a core can keep growing at the pressure max- 
imum, the core is required to reach ~ lOM^ within sev- 
eral Myrs for the formation of Saturn. For > 10 km, 
an MMSN disk {xs = = 1) cannot produce a core 



with 10 Earth masses in several million years (see Fig.|31). 
Fragment production from large planetesimals is ineffec- 
tive because of large and occasional collisions among 
them, which delays embryo growth via fragment accre- 
tion. A large solid surface density {xs > 1) that increases 
fragment production is necessary for planetesimals with 
larger rp to form a massive core in several million years. 
A core can reach 10 Earth masses in a disk with Xs > 3 
for ro — 10 km and with Xs ^ 10 for rg = 100 km (see 
Fig. m). Note that the growth of cores is almost inde- 
pendent of Xg, because fragments with rf contributing 
most to the growth are controlled by Stokes drag (e.g., 
iBeauge et al.lll994[ ). 




Time [year] 

Fig. 4. — Embryos can reach 10 Earth masses at 9AU in a 
massive disk with Xb = 3 for ro = 10 km and with Xb = 10 for 
ro = 100 km. The gas density of the disk is determined as Xg = 1 
(solid lines) and Xg = Xb (dotted lines). 

5. DISCUSSION 

In this paper, we have taken the size of initial plan- 
etesimals as a free parameter. The size depends on 
the formation mechanisms of planetesimals. If planetesi- 
mals are for med via self-gra vitat i onal i n stability in a thin 
dust layer (iGol dreich & War3 119731: iMichikoshi et al.l 
120121 : iTakeuchi fc Idal |2012[ ) or of dust accumulated in 
the structure of disk turbulence ([johanscn et al. 200^ 
ICuzzi et al.l 120081 ) ■ resulting planetesimals are 10 km or 
larger. On the other hand, coUisional coagulation of 
fluffy dust is possible to produce planetesimals over- 
comi ng an obstacle o f rapid radial drift at fstop 1 
( Oku zumi et al.l [20121 ): the resulting planetesimals may 
be smaller. In addition, bodies in a significant amount 
drift from the outer disk during the coUisional growth 
of fluffy dust and such drifting dust increases the solid 
surface de nsity form MMSN to about 4 x MMSN inside 
~ 10 AU (|Okuzumi et all [20121 ). Therefore, we should 
address the history of Saturn core formation including 
planetesimal formation. 

A planetary embryo born at the pressure maximum mi- 
grates due to interaction with the surrounding gas disk 
(type I) . The time required to go through half the width 
of the embryo's feeding zone due to type I migration is 
longer than the growth time of the embryo unless the 
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embryo exceeds OAM^ (jTanaka et al.ll2002D . The type 
I migration rate of an embryo with mass M at distance 
a is given by da/dt = 2fiMY,ga^^l^/ MqC^ , where /i is 
a factor that is determined hy p = — dlnEg/dhra and 
q — — dlnc/dlna. If horseshoe torques are fuhy satu- 
rated, /i = —1.45 — 0.6p — 1.8g in a locahy isothermal 
disk and /t = —1.04 — 0-77 p — Q.llq in an adiabatic disk 
(jPaardekooper et aLll20lil ). In the outer disk where we 
assume p = 1.5 and q = 0.25, type I migration is inward 
due to a negative vakie of /i, while outward migration 
occurs around the gap opened by Jupiter where the ra- 
dial Sg slope has a large positive value. Therefore type 
I migration stalls at a certain distance under the density 
profile shown in Fig. [TJ the location of zero migration is 
at 8.6 AU in the locally isothermal disk and at 9.1 AU in 
the adiabatic disk. Since the location is near the pressure 
maximum, such a large embryo, which has a wide feed- 
ing zone, can grow via the accretion of bodies accumu- 
lating at around the pressure maximum. Note that if the 
torques are fully unsaturated, /i = — 0.85— p-|-1.8(7 in the 
locally isothermal disk and /i = —0.61 — 2.33 p -I- 2.82(7 in 
the adiabatic disk ( Paardek ooper et al.ll2Q10f) . Then, /i 
is positive (outward migration) if q has a large value. In 
the optically thick regions (inner disks or massive disks) 
where viscous heating is more dominant than irradiation, 
q has a relatively large value. However, in the outer 
regions (~ 9AU) of moderate disks with ^ 1 that 
we are concerned with, it is likely that the disk temper- 
ature is determined by irradiation rather than viscous 
heating and that q has a relatively small value (~ 0.25). 
Then, for p = 1-1.5, type I migration is inward (/i < 0) 
and the argument of the migration trap is still appli- 
cable. Furthermore, the torque saturation is likely in 
the disk of interest. The saturation is prevented only if 
the horseshoe liberation time is comparable to the dif- 
fusion/ thennBJ_djfh]S^^ across the horseshoe 
width ([Paardekooper et al.ll201lD and if an ein bryo has a 
low enough eccentricitv (jBitsch fc Klevll2010[ ). The em- 
bryos resulting form our simulations have eccentricities 
larger than 0.01. Therefore the condition for type I mi- 
gration in the outer disk where Saturn's core forms is 
most likely to result in inward migration, except in the 
vicinity of the gap. Consequently, type I migration does 
not inhibit the fast formation of Saturn's core but rather 
strengthens it. 

Once a core exceeds a critical core mass, the core has 
no longer a static atmosphere and rapid gas accretion 
starts to form a gas giant. Interior modeling of Jupiter 
and Saturn suggests that Jupiter's core mass is smaller 
than lOMm, while Saturn's core is larger than lOM^ 
(IGuillotll200a . Since heating by the accretion of small 
bodies onto a core stabilises the atmosphere, the criti- 
cal core mass beco mes large for a fast growth core (e.g., 
Ilkoma et aLll2000D . Traditional scenarios, in which the 
core growth is 5-10 times slower for Saturn than for 
Jupiter, cannot explain the estimate of core masses of 
Jupiter and Saturn. However, Saturn core grows rapidly 
at the pressure maximum caused by Jupiter's gap open- 
ing. Since the growth is possible to be faster for Saturn 
than for Jupiter, the resulting core masses of Jupiter and 
Saturn may be consistent with those derived form the in- 
terior modeling. 

Since the core of Saturn is estimated to be larger than 



10 Me (|Guillotl 120051 ) ■ Saturn should be formed via core 
accretion. We have shown that if Jupiter has already 
been completed and opened a gap in the disk, the massive 
core of Saturn has rapidly formed from kilometer sized 
planetesimals in an MMSN disk or from larger planetes- 
imals in a disk with a larger amount of planetesimals. 
On the other hand, since small initial planetesimals are 
quickly removed by combination of fragmentation with 
gas drag, the core of Jupiter is likely to have formed from 
large initial planetesimals wi th ~ 100 km at a round 5 AU 
in a 10 X MMSN disk (Koba vashi et all 120 111) . Jupiter is 
possible to form earlier than Saturn core formation, if ei- 
ther (1) the solar nebula is as massive as self-gravitational 
instability occurs, (2) pressure maxima originally exist in 
the solar nebula, or (3) planetesimals of size 100 km are 
born around 5AU in a massive disk. All the cases are 
consistent with the sequential formation of Saturn. 

After the gas accretion, Saturn might open up another 
gap in the disk, resulting in the emergence of an addi- 
tional pressure maximum. Neptune and Uranus mainly 
consist of ice with more than lOM^ covered with tenuous 
atmospheres. Near the pressure maximum, such massive 
icy planets can be formed. Fragments produced in the 
outer disk accumulate within about the width of the disk 
scale height near the pressure maximum but the feeding 
zone of such a massive planet is wider than the scale 
height; hence only one planet can be produced. Mean- 
while, the orbital eccentricity distribution of Kuiper belt 
objects in the 3:2 mean-motion resonance with Neptune 
suggests that Neptune migra ted outwa,rd fro m around 
20 AU to the current position ()Malhotralll995[) . Neptune 
was therefore expected to have formed at the pressure 
maximum and to have then migrated outward due to in- 
teraction with the outer planetesimal disk (e.g.. llda et al.l 
2000; Kirsh ct al. 2009). Note that torques from inter- 
action with gas that cause type I migration suppress 
the planetesimal-driven migration. Since Saturn should 
form in gas with a high density, high gas-to-solid sur- 
face densities render type I migration stronger compared 
to the planetesimal-driven migration. Therefore, the 
planetesimal-driven migration is negligible for the Sat- 
urn formation. However, Neptune plausibly migrated 
outward through interactions with planetesimals because 
of a low gas density. Uranus might then have formed at 
around the pressure maximum caused by Saturn, where 
fragments produced from remnant planetesimals that in- 
duce the Neptune migration accumulate. Therefore, we 
should also address Uranus and Neptune formation after 
Saturn's formation. 

6. CONCLUSION 

We have investigated the core formation of Saturn at 
the pressure maximum caused by Jupiter's gap opening 
in the solar nebula. Small particles feel strong gas drag 
and drift radially. The loss of bodies due to the drift 
stalls emb ryo growth before embryos r each the critical 
core mass (jKobavashi et al.|[20Tol[20TTI ). Although bod- 
ies stirred by a large embryo can drift from the pressure 
maximum, the drift halts near the pressure maximum 
due to positive radial slope of gas pressure; hence an em- 
bryo at the pressure maximum effectively grows without 
loss of surrounding bodies (see Section [S]). Starting from 
monodisperse planetesimals, planetary embryos are gen- 
erated via coUisional evolution and stir remnant planetes- 
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imals, resulting in fragmentation of planetesimals. Frag- 
ments collide with each other and they accumulate at a 
radius rf, given by Equation The random velocity 

of bodies with rf is well damped by gas drag; thereby em- 
bryos rapidly accrete such bodies. Fragments produced 
in the outer disk move to the pressure maximum, which 
contribute to further embryo growth. Due to the rapid 
accretion, a core as massive as 10 Earth masses forms in 
several million years, for kilometer-sized initial planetes- 
imals in an MMSN disk, while the core formation needs 
a disk with 3 times larger solid surface density for 10 km 
planetesimals and with 10 times solid surface density for 
100 km planetesimals. The growth is almost indepen- 
dent of gas density. Since the rapid formation in a disk 
with moderate mass is consistent with insignificant type 
II migration of Jupiter and larger core of Saturn than 



that of Jupiter, Saturn's core may have formed in the 
pressure maximum after Jupiter opened up a gap in the 
solar nebula. 
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